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as  an  elastio-plastic  fracture  criterion,  the  report  suggests  a new  method  for 

its  laboratory  determination.  Advantages  of  the  new  method  include: 

unambiguous  continuity  with  K^;  full  compatibility  with  the  existing  standard 

c 


for  COD  testing;  and  an  ability  to  deal  with  slow  crack  growth  and  the  effects 
of  lateral  constraint.  The  report  concludes  with  an  outline  of  future  work 
and  detailed  instructions  for  the  determination  of  in  a 3pt  bend  specimen 

with  a crack  length  to  width  ratio  of  between  0.25  and  0.6.  . 
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NOMENCLATURE 

A in-plane  area  of  a body  (used  in  formulation  of  J) 

crack  length 
material  thickness 

constraint  factor  used  to  describe  effective  elevation  of  yield  stress 


a 

B 


in  3pt  bending 
E Youngs  modulus 

E^  E/(l-«^)  plane  strain,  > E plane  stress) 

f (8),  g. .(®)»  hi-sCe)*  dimensionless  functions  used  in  describing  crack  tip 
ij  ij 


P 

G 


h(-) 
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stress-strain  distributions 
ratio  of  load  point  displacement  q,  to  mouth  opening  displacement  V 
available  energy  release  rate  in  a unit  thickness,  linear  elastic, 
cracked  body 


polynomial  used  to  describe  mouth  opening  compliance 


J 


J 
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k 


n 

P 

Q 

S’ 

q 

r 

r,e 

s 

s 

T 

U 

^ crack 


J integral,  or  rate  of  change  in  absorbed  energy  per  unit  thickness 
with  crack  length  for  specimens  loaded  monotonically  to  a given  applied 
displacement 

value  of  J at  unstable  fracture 

elastic  and  plastic  components  of  J 

linear  elastic  stress  intensity  factor 
critical  value  of  K in  plane  strain 

value  of  K at  final  load 

plastic  stress  intensity  factor 
work  hardening  exponent 
potential  energy 
load 

limit  load  or  net  section  yielding  load 
final  load 

load  point  displacement 

constant  used  to  describe  centre  of  rotation  during  plastic  deformation 
crack  tip  polar  co-ordinate  system 
distance  along  perimeter  of  cracked  body 

distance  along  a contour  path  (used  in  formulation  of  J integral) 
external  boundary  forces 

area  under  the  load/load  point  displacement  curve,  or  total  absorbed 
energy 

component  of  U due  to  the  presence  of  the  crack 


V 


u , u 

e p 

elastic  and  plastic  components  of  U 

V 

mouth  opening  displacement 

V 

e 

displacement  of  knife  edges  positioned  at  crack  mouth 

w 

width  of  specimen 

w 

strain  energy  density 

x,y 

cartesian  co-ordinate  system 

Y 

geometry  factor  used  in  calculating  K 

z 

height  of  knife  edges  above  specimen  surface 

Greek  Symbols 

a 

constant  used  in  work  hardening  relationship 

P 

geometry  dependent  constant  used  in  slow  crack 

growth  correction 

r 

crack  tip  integration  path  (used  in  formulation  of  J) 

5 

crack  opening  displacement  (COD) 

6c 

critical  value  of  6 (at  unstable  fracture) 

6 

P 

plastic  component  of  6 

< 

strain 

€P 

plastic  strain 

0 

angular  co-ordinate 

geometry  dependent  constants  used  in  linking  J 

to  absorbed  energy 

a 

stress 

a 

y 

yield  stress 

» 

poisson's  ratio 

Abbreviations 

COD 

crack  opening  displacement 

CTS 

compact  tension  specimen 

LEFK 

linear  elastic  fracture  mechanics 

THE  J INTEGRAL.  A REVIEW  OF  THEORY  AND  A SUGGESTED 

METHOD  FOR  J DSTERiTN ATT  ON 
c 

by 

JPG  Sumpter 

INTRODUCTION 

1.  The  objective  of  fracture  mechanics  analysis  is  to  predict  the  behaviour 

of  defects  in  a large  complex  structure  from  the  results  of  a small,  easily 

performed,  laboratory  test.  Most  attempts  to  do  this  are  based  on  the  concept 

of  a fracture  criterion,  the  critical  value  of  which  is  postulated  to  be  the 

same  in  both  test  piece  and  structure.  The  two  most  widely  used  fracture-safe 

design  procedures  are  those  based  on  linear  elastic  fracture  mechanics  (LEFM), 

fracture  criterion  K..  , and  on  the  crack  opening  displacement  (COD)  concept, 

1c 

fracture  criterion  6 . 

c 

2.  Strict  application  of  LEFM  requires  that  crack  tip  plasticity  be 
negligible  in  both  test  piece  and  structure.  This  requirement  is  met  in 
practice  by  stipulating  that  (l)# 


a,  B,  W-a  » 


Ul 


where  a is  crack  length,  B is  material  thickness,  W is  specimen  width  and  a 


is  the  material  yield  strength.  The  ductility  of  most  engineering  materials 
is  such  as  to  make  this  condition  an  extremely  restrictive  one,  and  it  is 
seldom  possible  to  obtain  a valid  K^c  at  the  thickness  and  temperature  of 


interest.  Most  fracture  analyses  must  consequently  be  based  on  a parameter 
which  can  be  measured  in  a small  test  piece  irrespective  of  the  extent  of 
plasticity  at  fracture. 

3.  The  COD  approach  fulfils  this  requirement.  The  basic  assumption  of  the 

method  is  that  a crack  deforms  with  a well  defined  opening  at  its  tip  and  that 

fracture  occurs  when  this  discrete  opening  reaches  a critical  value  £>  • A 

c 

standard  method  for  determination  of  6 is  laid  down  in  (2).  The  test  piece 

c 

is  conventionally  a 3pt  bend  specimen,  span  to  width  ratio,  S/W  - 4,  crack 


length  to  width  ratio,  0.25<a/w<0.6,  and  thickness  to  width  ratio  B/W  = 1.0  or 
B/W  . 0.5.  COD  is  inferred  from  a plot  of  load  against  displacement,  the 
latter  being  monitored  by  a clip  gauge  mounted  across  the  mouth  of  the  crack. 

4.  Although  determination  of  6^  from  a laboratory  test  piece  is  reasonably 

straightforward,  the  use  of  COD  in  fracture-safe  design  is  less  satisfactorily 
established.  The  main  difficulties  arise  from  an  inability  to  provide  accurate 


*(  ) - References  on  page  14-15 


COD  solutions  when  geometry  i3  complex  and  stress  gradient  significant;  and  from 
the  uncertain  relationship  between  6 and  X for  near  elastic  behaviour.  Suggested 
procedures  for  the  use  of  COD  in  critical  defect  size  analysis  (3,  4)  overcome 
these  difficulties  by  the  use  of  a universal  design  curve.  This  inevitably 
contains  a degree  of  conservatism  in  some  cases  and  there  are  advantages  to 
be  gained  from  using  a fracture  criterion  which  can  provide  explicit  solutions 
in  terms  of  crack  size  and  stress  distribution  at  any  structural  location.  The 
J integral  fulfils  this  requirement  and  has  the  added  advantage  of  providing 
unambiguous  compatibility  with  X for  near  elastic  behaviour  (paragraph  8).  The 
use  of  J for  fracture  design  is  further  discussed  in  (p,  6,  7) • 

5.  No  generally  agreed  procedure  yet  exists  for  determination  of  J although 

c 

there  is  considerable  interest  in  its  use,  especially  in  the  US,  where  the  E24 

committee  is  formulating  a draft  standard.  This  report  describes  a method  for 

obtaining  J using  the  instrumentation  and  specimen  geometry  already  standardised 
c 

for  COD  testing.  COD  practice  is  followed  further  in  defining  J as  the  point 

c 

of  unstable  fracture  (which  may,  or  may  not,  be  preceded  by  stable  crack  growth) 
in  a test  piece  having  the  same  thickness  as  the  structure  of  interst. 

THEORETICAL  BASIS  FOR  THS  USS  OF  J A3  A FRACTURE  CHARACTERISING  PARAMETER 

6.  It  is  reasonable  to  assume  that  fracture  initiation  from  a sharp  crack  is 
governed  by  the  very  high  stresses  and  strains  which  exist  ahead  of  that  crack. 

An  ability  to  quantify  this  stress-strain  field  in  terms  of  overall  geometry  and 
loading  should  consequently  lead  to  an  ability  to  predict  the  conditions  under 
which  fracture  will  occur.  It  is  well  known,  that  having  made  the  assumption 

of  perfectly  linear  elastic  behaviour,  it  is  possible,  on  a continuum  mechanics 
scale  at  least,  to  provide  just  such  a complete  description  of  the  near  crack 
tip  region. 


K 


[2  ] 


ffij’  €ij  V2^  Zij 

where  a..,  «. . are  the  cartesian  co-ordinates  of  stress  and  strain  at  the  crack 

1 J 1 J 

tip;  r,  8 is  a polar  co-ordinate  system  centred  at  the  crack  tip;  and  K is  the 
linear  elastic  stress  intensity  factor.  Since  f ^ - (© ) is  unique  for  r->  0 (8), 
only  the  magnitude,  X,  of  the  1/Vr  singularity  is  dependent  on  geometry  and 
loading  conditions.  X consequently  provides  a single  parameter  description  of 
the  crack  tip  environment,  and  its  critical  value  may  be  expected  to  provide  a 
geometry  invariant  fracture  criterion,  as  long  as  the  assumption  of  linear 
elastic  conditions  at  the  crack  tip  is  a valid  one.  As  noted  in  paragraph  2, 

K has  in  fact  been  shown  to  have  a geometry  invariant  critical  value,  K1(j, 
provided  the  size  constraints  of  [1]  are  met.  Outside  these  constraints  the 
size  of  the  cr^ck  tip  plastic  zone  is  such  as  to  invalidate  [2]  and  the  critical 
value  of  K i3  no  longer  a material  property. 
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7.  In  principle,  it  should  he  possible  to  extend  the  continuum  mechanics 
description  of  the  crack  tip  region  to  account  for  plastic  deformation;  but, 
in  practice,  elastic-plastic  solutions  are  difficult  to  obtain  and  complete 
solutions  are  available  only  for  some  very  specialised  cases.  For  the 
important  case  of  a crack  under  a remotely  applied  tensile  stress  field  only  an 
approximate  solution  exists  (9,  10,  ll).  This  gives  the  near  tip  (r^  o) 
stress  distribution  well  inside  the  plastic  zone  of  a work  hardening  material 
following  a deformation  plasticity  flow  rule  with  no  unloading  (essentially 
non-linear  elasticity  in  that  stress  and  strain  have  a unique  relationship) . 

No  limitation  is  placed  on  the  size  of  the  plastic  zone,  but  finite 
boundary  conditions  are  not  explicitly  accounted  for.  The  important  point 
to  emerge  is  that,  once  again,  as  for  the  linear  elastic  case,  there  is  a 
unique  distribution  of  singular  terms  at  the  crack  tip.  For  a material 
following  a stress-strain  behaviour  of  the  form 
P 


r /■  \n 

€ J 2L1 


[3] 


the  near  tip  stress-strain  distribution  is  of  the  form  (10,  ll) 
-1 


-n 


L* 


* n | x 
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(n+l) 


V8) 


where  c is  plastic  strain,  n is  the  work  hardening  exponent,  k is  a plastic 
stress  intensity  factor,  and  g^  (0)  and  h^  (©)  are  dimensionless  functions 

detailed  in  (10,  ll).  Although  correlations  for  k have  been  provided  in  some 
configurations  (12)  using  specially  formulated  finite  elements,  it  is  more 
convenient  to  have  the  strength  of  the  dominant  singularity  expressed  in  terms 
of  a parameter  more  directly  obtainable  from  the  overall  load  deflection 
behaviour  of  the  cracked  body.  The  J integral  provides  a particularly  con- 
venient method  of  doing  this . 

t).  For  a body  of  unit  thickness,  composed  of  a linear  or  non-linear  elastic 
material,  the  available  energy  release  rate  or  crack  extension  force  G is 
defined 


-dP 

da 
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ra 


f T.  du,  dS  - 

f dwdA~l 

Li^  1 j 

'a  J 

[5] 


where  P is  the  potential  energy,  T^  and  u,  are  the  forces  and  displacements  on 
the  body  perimeter  S,  A is  the  area,  and  w is  the  strain  energy  density. 
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Suppose  the  body  is  described  by  a co-ordinate  system  x,y 5 by  a transformation 
of  co-ordinates,  most  easily  seen  for  the  steady  state  case  where 


-dP  _ dP  I 

da  = dx  I y,a  [6] 

together  with  application  of  Green's  theorem  to  the  term  j ,dwdA,  the  HHS  of  [ p] 

can  be  expressed  as  a path  independent  line  integral,  designated  J by  Rice  (l 3) * 


G 


. T.  du.  ds 

jr(wdy " '“S  }-J 


where  s is  the  distance  along  any  arbitrary  contour  path,  r,  surrounding  the 
crack  tip.  It  follows,  that  for  linear  elasticity,  the  identity  between  G and 
K (8)  holds  also  between  J and  K. 


K2 

J - % [8] 

E 

where  E^  « E plane  stress 
E 

» -■  ■ plane  strain 

(l-»2)  ‘ 

u is  Poisson's  ratio. 

9.  The  potential  energy  decrease  which  occurs  during  crack  extension,  A a,  may 
be  conveniently  visualised  as  the  area  between  the  overall  load  deflection 
curves  of  a body  containing  a crack  length  a,  and  that  of  a body  containing  a 
crack  length  (a+Aa).  It  can  be  shown  that  this  area  is  the  same  irrespective 
of  whether  load  or  displacement  are  held  constant  during  crack  extension  ( 1 3) . 
This  being  so,  [5]  can  be  simplified  by  taking  the  energy  change  under  fixed 
boundary  displacements.  The  change  in  external  work/  Tdu. ds  is  then  zero,  and 

f 5 -1 

only  the  stored  energy,  U = i^wdA,  changes  (figure  l) . For  the  general  case  of 


a body  thickness  B 


where  q.  is  the  load  point  displacement.  For  a plastic  material,  where  U is 
dissipated  in  irreversible  plastic  work,  J can  no  longer  be  identified  with  the 
elastic  energy  release  rate  G (14) , but  it  will  still  be  the  rate  of  decrease 
of  area  under  the  load/load  point  deflection  curve  with  crack  length.  This 
makes  it  a relatively  easy  parameter  to  determine  experimentally  (paragraph  15). 
10.  The  J integral  with  w and  u representing  total  (elastic  plus  plastic) 
strain  energy  density  and  displacement  will  also  retain  path  independence  for  a 
material  following  a deformation  plasticity  flow  rule.  This  allows  it  to  be 
regarded  as  a measure  of  the  crack  tip  singularity  revealed  in  (9,  10,  ll) 
Evaluation  of  J as  n*  0 may  be  achieved  by  substitution  of  [4]  into  [7],  It 
follows  that 
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[10] 
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where  I is  a dimensionless  ccnstant  (10,  ll).  [10]  provides  the  basis  for  an 

elastic-plastic  fracture  mechanics  approach  analogous  in  form  to  IEFK.  The 

hypothesis  is:  that  there  is  a unique  stress-strain  distribution  at  the  crack 

tip;  that  J characterises  this  stress-strain  distribution;  and  that  J might 

c 

consequently  be  expected  to  provide  a geometry  invariant  fracture  criterion. 

11.  It  should  be  emphasised  that  the  above  hypothesis  can  be  criticised  on 
a number  of  grounds.  Firstly,  the  derivation  of  [10]  is  based  on  deformation 
flow  rule  plasticity.  Real  materials  are  better  represented  by  an  incremental 
flow  rule  (where  strain  is  history  dependent,  not  uniquely  defined  by  the 
current  value  of  stress).  Differences  between  the  two  theories  are  only  large 
when  significant  rotation  of  the  principal  stress  system  or  unloading  occurs, 
but  this  cannot  be  ruled  out  near  the  crack  tip,  especially  if  slow  tearing 
precedes  final  fracture.  Secondly,  available  slip  line  field  solutions  for 
different  cracked  body  geometries  would  seem  to  indicate  that  a single  para- 
meter description  of  the  crack  tip  may  no  longer  be  possible  in  the  presence 
of  different  patterns  of  plastic  defoliation  (l5)»  It  has  been  argued  (l6,  17) 
that  work  hardening  will  overide  these  differences  and  create  a unique 
environment,  or  automonous  zone,  at  the  crack  tip.  Elastic-plastic  finite 
element  analysis,  in  so  far  as  it  can  be  expected  to  provide  accurate  answers 
in  the  near  tip  region,  does  not  seem  to  support  this  view  (l8,  19). 

12.  A number  of  questions  thus  surround  the  theoretical  basis  for  using  J as 
a crack  tip  characterising  parameter  and  ultimately,  its  usefulness  as  a 
fracture  design  tool  is  likely  to  be  decided  by  experimental  evidence.  Before 
this  can  be  done  there  is  an  obvious  need  to  fix  on  a generally  agreed  pro- 
cedure for  J testing.  There  is  moreover  a requirement  that,  in  order  to 

c 

encourage  its  widespread  use,  any  proposed  test  method  should  be  reasonably 

easy  to  perform  and  preferably  fully  compatible  with  existing  standards  for 

K and  6 testing  (l,  2).  The  remainder  of  this  report  is  devoted  to 
1.  c c 

suggesting  such  a method. 


DETERMINATION  OF  J 

c 

Use  of  Finite  Element  Based  Correlations 

13.  The  contour  integral  definition  of  J [7]  is  well  adapted  to  evaluation 
from  elastio-plastic  finite  element  analysis  (20)  and  a large  number  of  studies 
have  been  performed  on  test  piece  geometries.  Such  computations  provide 
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correlations  between  J and  apniied  displacement  which  sho  id  be  directly  stnlic- 
nble  to  fracture  toughness  testing.  Figure  2 shows  plane  strain  (2l)  and  plane 
stress  (22)  results  for  J against  mouth  opening  displacement  V in  the  3pt  bend 

geometry  as  an  example.  Unfortunately,  these  correlations  are  rather  inflexible 
since  each  one  refers  only  to  a specific  geometry  and  does  not  in  itself  give 
any  indication  of  how  the  value  of  J will  be  affected  by  minor  variations  away 
from  this  configuration.  Crack  length  in  particular  cannot  be  specified 
a-priori  in  standard  fracture  toughness  tests.  It  can  moreover  be  shown  that 
the  various  estimation  procedures  described  below  (paragraph  14)  all  show  close 
agreement  with  computed  results  for  J (23)*  This  together  with  their  greater 
flexibility  means  that  they  are  more  suitable  than  finite  element  based  correl- 
ations for  determination  of  J in  standard  test  piece  geometries.  1'urcerical 

c 

analysis  still  has  a major  role  to  play  in  the  analysis  of  non-standard  test 
pieces  such  as  those  containing  very  shallow  notches,  and  for  investigation  of 
unusual  structural  configurations  (23)« 

Correlations  Between  J and  Area  under  the  Load  Deflection  Curve 

14.  The  most  convenient  basis  for  estimation  of  J in  small  test  pieces  is  that 

provided  by  [ 9] , where  for  an  elastic-plastic  material  U is  the  total  absorbed 

energy,  or  area  under  the  load/load  point  deflection  curve.  Early  attempts  to 

establish  J (16,  24)  used  a series  of  specimens  with  different  initial  notch 
c 

depth  to  evaluate  J directly  from  [ 9] » the  possibility  of  obtaining  Jc  from  a 

single  specimen  load  deflection  curve  using  an  estimation  procedure  based  on 
elastic  compliance  and  slip  line  field  considerations  was  Inter  pointed  out  (25) t 
and  more  recently  (26)  attention  has  focussed  on  the  use  of  correlations  between 
J and  the  total  absorbed  energy.  The  first  of  these  (27),  for  specimens  in 
which  the  uncracked  ligament  is  subjected  to  bending,  took  the  form 

T 2 U crack  r , , n 


0 ■ (W--a)  T L^j 

where  ^cracj<;  denotes  the  component  of  energy  absorbtion  due  to  the  presence  of 

the  crack.  For  deeply  cracked  (a/w  > 0.6)  3pt  bend  and  compact  tension 

specimens,  where  the  energy  absorbtion  component  of  the  uncracked  body  is  small 
it  was  suggested  (26)  that  [11]  be  simplified  to 

where  U is  once  again  the  total  area  under  the  load/load  point  deflection  curveo 
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15.  A more  general  version  of  [ 11 J , which  does  not  involve  subtraction  of 
the  uncracked  body  energy  and  is  applicable  to  any  geometry  for  which  the 
elastic  compliance  and  limit  load  are  known,  can  be  written  in  the  form  (23) 


J = J + J 
e p 


T)  U T)  U 
e e op 
+ 


J and  U are  here  divided  into  elastic  and  plastic  components  (figure  3).  In 
contrast  with  [11],  U here  includes  the  energy  of  the  uncracked  body,  tj  an 

T)  are  geometry  dependent  constants.  T}  is  based  on  the  stress  intensity 
P 6 

factor,  which,  in  turn,  can  be  used  to  calculate  J from  [8],  and  U from 

© © 

load  point  compliance  relationships  (28).  Table  1 shows  t)  values  for  the 
3pt  bend  and  compact  tension  geometries  based  on  stre -,s  intensity  and  com- 
pliance relationships  given  in  (25)  and  (29)  respectively.  r|  may  be  cal- 

p 

culated  from  a knowledge  of  the  relationship  between  crack  length  and  limit 
load  For  example,  for  the  3P"t  bend  specimen,  provided  that  plasticity 

is  confined  to  the  uncracked  ligament,  the  net  section  yielding  (or  limit) 
load  is  a function  of  ligament  length  squared 

- Cn  ffy  5 (w"a)2  [14 

where  is  a constant.  For  a non-work  hardening  material 


. r dU  \ 

JP  - h{  tt) 

UP  - QLqp 


whence 


__  (W- a)  _ 

2C  a > ‘ q_ 

n v 5 t) 


p 


It  follows  that  ri  = 2 for  the  3P^  bend  geometry  at  all  crack  to  width  ratios 

greater  than  about  a/W  * 0.2.  Limit  load  for  the  CTS  is  less  well  defined, 
but  tj  might  be  expected  to  approach  2 for  deeply  cracked  specimens  where  the 

defoliation  is  predominantly  bending. 


TASLS  1 

VAIL'23  CF  tj  POh  [13] 


16.  For  the  CTS  (a/w  > C.6),  and  the  3pt  bend  geometry  (0.4  < a/w  < 0. 7)  where 
Tl . -T,  n - 2,  [13]  reduces  to  [12].  It  should  be  re-emphasised,  however,  that  the 

0 p 

values  of  n in  Table  1 do  not  involve  subtraction  of  an  uncracked  body  energy. 

'e 

There  will  therefore  be  a discrepancy  between  [13]  and  [ll]  when  the  uncracked 
body  energy  is  significant,  as  it  is  in  the  3pt  bend  geometry  (30?1  of  the  total 

elastic  at  a/w  . 0.5,  70%  at  a/w  = 0.3).  It  has  recently  been  suggested  (30)  that 
the  source  of  this  discrepancy  is  an  error  in  the  formulation  of  [11],  but,  at 
the  time  of  writing  the  question  has  still  to  be  completely  resolved.  The  draft 
US  standard  for  J testing  (31)  seems  likely  to  adopt  the  specimen  design  and 
procedure  suggested  in  (26),  namely  a deeply  notched,  (a/w  > 0.6)  CTo  -.nth  the 
mouth  of  the  notch  machined  away  to  allow  for  insertion  of  the  clip  gauge 
directly  under  the  loading  pins.  This  type  of  specimen  is  not  always  satis- 
factory, especially  when  testing  weldments  (32)  and  it  remains  to  be  seen 
whether  the  proposed  specimen  design  is  revised  in  the  light  of  (30). 

17.  Application  of  [13],  or  its  simplified  version  [12],  requires  measurement 
of  load  point  displacement.  This  is  not  a quantity  routinely  measured  in 
fracture  toughness  testing.  The  standard  procedure  (l,  2)  is  to  measure  crack 
mouth  opening  V.  As  noted  above,  it  is  possible  for  the  CTS  geometry,  to 
relocate  the  clip  gauge  directly  under  the  loading  pins,  thereby  measuring  load 
point  displacement  directly.  This  option  is  not  available  for  the  3pt  bend 

geometry.  Rather  than  attempt  a direct  measurement  of  load  point  displacement 
in  this  specimen  it  is  probably  preferable  to  convert  V to  load  point  displace- 
ment q,  by  assuming  a centre  of  rotation  a fixed  distance  r(W-a)  below  the 
crack  tip.  So  that,  before  application  of  [l3]»  the  area  under  a load/mouth 
opening  deflection  curve  should  be  factored  by  F,  where 


.VS* 


It  is  well  known  from  studies  in  support  of  the  COD  concept  (33)  that  r is 
both  load  and  geometry  dependent.  After  limit  load,  however,  r does  achieve 
a fixed  value  for  increments  in  q and  V in  a given  geometry.  Finite  element 
results  (22)  for  the  3pt  bend  geometry  in  plane  strain  indicate  that  this 

fixed  value  is  0.45  at  a/tt  * 0.3,  and  0.4  at  a/V  = 0.5.  The  plastic  com- 
ponent of  F can  be  calculated  from  [19]  using  these  values.  Modifying  [19] 
for  use  with  areas  under  a load/mouth  opening  displacement  curve  by  employing 
separate  values  of  F for  elastic  and  plastic  behaviour  results  in 

F n U F n U 
T e 'e  ve  p 'p  vp 

J - (tf-aiif  + -tn-fjr-  [20] 

where  U and  U are  the  elastic  and  clastic  areas  under  the  load/mouth 
ve  vp 

opening  deflection  curve.  Relevant  values  for  the  3pt  bend  specimen 

a/w  a 0.3  and  0.5,  and  the  CTS  a/w  . 0.5  are  summarised  in  Table  2.  F for 

the  3pt  bend  geometry  comes  from  boundary  collocation  results  for  elastic 

mouth  opening  given  in  (34);  F for  the  CTS  comes  from  (35);  2nd  F for  the 

0 P 

CTS  is  calculated  on  the  basis  r * 0.4  (33). 

TABLE  2 

GEOMETRY  FACTORS  FOR  USE  WITH  T 201 


Specimen 

3pt  bend  a/W  - 0.3 

3pt  bend  a/w  = 0.5 

CTS  a/w  =0.5 

F 

e 

2.45 

I.65 

0.70 

F T) 
e ‘e 

4.20 

3.25 

1.60 

F 

P 

1.65 

1.45 

0.75 

F T) 

PP 

3.25 

2.85 

1.50 

A New  Method  for  J Determination,  based  on  Load,  and  Clip  Gauge  Plastic 
Displacement 

18.  Although  [20]  provides  a reasonably  satisfactory  method  of  determining 
J it  is  worth  questioning  whether  measurement  of  absorbed  energy  really 
provides  the  best  basis  for  a new  standard  test  method.  The  two  fracture 
toughness  tests  at  present  in  most  widespread  use  are  based  on  load  (K^) 

and  clip  gauge  measured  mouth  opening  (6  ).  There  is  much  to  be  said  for 

c 

providing  as  close  a compatibility  as  possible  with  these  established  test 
procedures.  At  the  linear  elastic  end  of  the  range  particularly,  it  is 
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essential  that  there  should  be  continuity  between  Jc  and  K^.  Discrepancies 

between  experimental  and  theoretically  predicted  elastic  compliances  can  easily 
upset  this  agreement  if  u ^ is  determined  in  the  usual  manner  from  load  whilst 

J is  determined  from  area  under  the  curve.  There  seems  to  be  a strong  case 
c 

for  forcing  compatibility  between  K and  the  elastic  component  of  J by  use  of 
[3]  . Suppose  the  load  at  failure  (figure  3)  is  Then 


[21] 


[22] 


where  Y is  a function  of  a/w  (l). 

19.  Since  Jg  is  to  be  evaluated  from  load  in  accordance  with  (l)  it  seems  con- 
sistent to  follow  GOD  testing  practice  (2)  and  determine  Jp  from  mouth  opening 
displacement.  This  can  be  done  by  a straightforward  substitution  of  a from 


[19]  in  [17] 
J 


2 C a F V 

n y S p p 


[23] 


Evaluation  of  C could  be  based  on  slip  line  field  theory.  For  3pt  bending  the 
n 

value  of  this  constant  is  about  1.5  (36,  37) > but  a better  alternative,  which 

does  not  involve  the  assumption  that  the  analytical  and  experimental  limit  loads 

will  coincide  is  to  substitute  for  C from  [14] 

n 

V5 


Bo  (V-a)2 

J 


whence 


[25] 


20.  Strictly  speaking,  [25]  applies  only  to  a non-work  hardening  material, 
where  is  both  the  load  at  net  section  yielding  and  the  maximum  load  reached. 

In  this  paper  is  used  to  signify  the  net  section  yielding  (not  the  maximum) 


load  in  a work  hardening  material.  In  most  geometries  and  materials  this  is 
denoted  by  an  easily  definable  inflexion  in  the  experimental  load  deflection 
curve.  [13]  indicates  that  J at  a given  deflection  will  be  higher  in  a work 
hardening,  than  in  a non-work  hardening  material,  whilst  [25]>  unless  it 
includes  a slight  rise  in  net  section  yielding  load,  will  not  give  this  trend. 
Evidence  on  the  effect  of  work  hardening  on  J is  not  conclusive.  The  J estim- 
ation procedures  for  cracked  test  pieces  in  tension  presented  in  (27)  all  pre- 
dict an  increase  in  J with  work  hardening,  but  finite  element  analysis  on  this 
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geometry  (23)  does  not  seem  to  exhibit  this  trend.  Finite  element  results 
for  3pt  bending  show  a similarly  small  effect  (22).  It  would,  of  course,  be 

fairly  simple  to  replace  in  [25]  by  (9^  + ^)/2,  or  even  9^,  but,  at  present, 

it  seems  more  sensible  to  follow  the  usual  lower  bound  philosophy  of  fracture 
toughness  testing  and  retain  [25]. 

21.  It  is  noted  in  (26)  that  [12]  over-estimates  J when  slow  crack  growth 
occurs.  This  point  is  illustrated  in  figure  4 in  which  a^  denotes  the  crack 

length  at  the  start  of  the  test,  and  a^  the  crack  length  after  slow  growth  Aa. 

In  so  far  as  J can  be  said  to  have  any  significance  to  crack  tip  conditions 
after  slow  growth,  the  absorbed  energy  U in  [12]  should  properly  be  that  appro- 
priate to  a crack  length  a^  monotonically  loaded  to  a given  deformation  (shaded 

area  in  figure  4),  not  the  area  under  the  experimental  load  deflection  curve, 
which  represents  a crack  length  a.^  growing  to  length  a^  during  the  course  of 

the  test.  From  [23],  the  value  of  J after  slow  growth  is 
(W-a  ) 

J = 2C  a -y..-  F V [26] 

p n y 5 p p 

The  constant  in  [26]  can  be  evaluated  as  before  from  [24] , but,  assuming  no 
slow  crack  growth  before  limit  load,  the  relevant  crack  length  is  a^  not  a^ 

QLS 

C 2 [27] 

V*"i> 

29L(K-aJ  F V 

J - — [3] 

P B(W-a.)2 

If  there  is  no  slow  growth  [26]  reduces  to  [25 ]•  If  slow  growth  does  take  place 
[28]  does  not  over-estimate  J in  the  way  that  [ 1 2]  does.  In  small  specimens, 
where  slow  growth  causes  a significant  change  in  elastic  compliance,  the 
measured  value  of  should  be  adjusted  appropriately.  This  procedure  is 

described  more  fully  in  the  Appendix.  F ahd  J should,  of  course,  be  evaluated 

P 6 

on  the  basis  of  a^.  not  a^. 

22.  As  shown  in  figure  2,  there  is,  for  3pt  bending,  a significant  differ- 
ence between  the  plane  stress  and  plane  strain  values  of  J.  This  difference  is 
less  pronounced  in  some  other  geometries,  notably  centre  cracked  tension  (23), 
but  in  general,  use  of  a finite  element  correlation,  or  of  an  estimation  pro- 
cedure of  the  type  suggested  in  (25)*  requires  a knowledge  of  whether  the  test 


piece  is  in  plane  stress  or,  pl?.ne  strain,  or  some  int  mediate  condition.  The 
use  of  [13]  or  [ 28 J circumvents  this  problem  since  the  lower  value  of  J in  plane 
stress  automatically  follows  from  the  lower  experimental  limit  load.  The  only 
assumption  required  is  that  the  limit  load  varies  as  the  square  of  the  ligament 
length  in  bending  (or  directly  with  the  ligament  length  for  pure  tension).  In 
passing,  it  should  be  noted  that  finite  element  analysis  (23)  indicates  that 
an  equal  value  of  J in  plane  stress  or  plane  strain  by  no  means  ensures 
identical  crack  tip  conditions. 

RELATIONSHIP  BETWEEN  J AND  GOD  IN  3pt  BENDING 


23«  There  are  two  recommended  methods  for  evaluation  of  COD  from  mouth 
opening  displacement  (2).  The  first  assumes  a fixed  centre  of  rotation  a 
distance  r(w-a)  below  the  crack  tip  throughout  the  loading  range 
r(w-a) 


r_r(w-a)  1V 

[a  + r(w-a;J 


[29] 


r is  put  equal  to  0.33  on  the  basis  of  experimental  investigations  (33)  to  give 
lower  bound  values  of  6.  The  second  method  due  to  Wells  (3B)  uses  an  estimation 
procedure  very  similar  to  that  proposed  for  J in  ( 25 ) • For  elastic  behaviour 
a relationship  6 = G/2.1 a is  assumed,  G being  linked  to  V through  elastic  com- 

y 

pliance;  at  an  assumed  limit  load,  derived  from  slip  line  field  theory,  the 
relationship  is  changed  to  [29],  but  with  r chosen  as  0.45,  again  on  the  basis 
of  slip  line  field  theory. 

24.  There  is  an  obvious  similarity  between  COD  defined  by  [29]  and  defined 

by  any  of  the  mouth  opening  displacement  related  equations  suggested  above,  in 
that  both  assume  a centre  of  rotation  a fixed  distance  r(w-a)  below  the  crack 
tip.  Considering  for  the  moment  the  plastic  component  of  COD,  5^,  [25]  and  [29] 

together  give 

J 2^W 

t*  - — P [30] 

p rB  (w-a) 

or,  assuming  a specific  value  of  r = 0.4  and  Cn  * 1.5,  [23]  and  [29]  give 

- Vp  t31) 

This  relationship  has  been  widely  noted  in  the  literature  from  finite  element 
and  other  analysis  on  the  3pt  bend  geometry  (21,  39>  40).  The  existence  of 


a similar  relationship  between  J and  6,  for  plane  strain  small  scale  yielding 
has  been  demonstrated  (40,  4l)  and  has  also  been  suggested  on  a more  intuitive 
basis  by  Wells  (42).  The  assumption  of  6 - G/2.1 o in  the  Wells  COD  formula 

•r 

follows  from  this  latter  reference  and  means  that  near  compatibility  should 
exist  between  total  values  of  J and  6 defined  in  this  way.  There  are,  however, 
a number  of  factors,  such  as  the  non-agreement  of  the  experimental  limit  load 
with  the  Wells  predicted  limit  load  which  may  upset  this  agreement  in  practice. 
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The  relationship  between  total  values  of  J and  COD,  defined  with  r = l/3  in 
the  elastic  and  plastic  regimes,  will  in  general  be  load  dependent,  varying  from 
J > 0.5 a 6 just  outside  the  ASTM  limit,  to  J > 'h  6 well  after  net  section 

y y 

yield  (23). 

SUMMARY  AND  DISCUSSION 

25.  The  proposed  method  for  determination  of  for  a 3?t  bend  specimen  is 

summarised  in  the  Appendix.  Few  of  the  ideas  used  to  derive  this  procedure 
are  completely  new  and  can  be  found  in  various  forms  in  previous  literature 
on  the  subject  (25,  27).  It  is  considered,  however,  that  the  method  proposed 
here  is  sin  especially  appealing  one,  both  in  its  simplicity  for  general 
laboratory  use  and  because  of  its  compatibility  and  similarity  with  existing 
fracture  toughness  testing  stsindards  (l,  2).  This  similarity  is  an  obvious 
advantage  in  gaining  acceptance  for  a new  testing  method.  Especially  in  the 
UK,  where  experience  and  confidence  has  been  gained  in  the  use  of  COD,  there 
is  an  understandable  reluctance  to  make  radical  changes  in  procedure.  As  shown 
here,  may  be  obtained  as  easily  as  6Q,  using  identical  test  pieces  and 

instrumentation.  The  two  major  reasons  for  favouring  the  use  of  J over  6 are 
its  unambiguous  continuity  with  K for  near  elastic  behaviour,  and  its  greater 
flexibility  and  ease  of  definition  in  fracture-safe  design  procedures  (5,  6,  7). 
FUTURE  WORK 

26.  A fairly  extensive  experimental  program  is  needed  to  fully  evaluate  the 

proposed  method.  Such  a program  would  need  to  cover  the  effect  of  specimen 

size  and  geometry  on  J , and  slow  crack  growth  prior  to  J . The  main  point 

c c 

of  comparison  would  be  with  the  US  proposed  method  ( 31 ) based  on  [12],  but 
comparisons  with  COD  would  also  be  made  to  establish  a link  with  UK  practice. 

27.  One  problem  which  has  already  come  to  light  is  the  very  early  onset  of 
slow  crack  growth  in  weld  metals.  In  large  specimens  (B  . 50  mm),  slow  tearing 
apparently  begins  well  before  the  limit  load  (43) » and  the  resulting  non- 
linearity obscures  the  well  defined  inflection  which  would  otherwise  occur  in 
the  load  deflection  curve  at  Q^.  This  makes  selection  of  an  appropriate  limit 

load  value  for  use  in  [26]  a difficult  procedure.  It  is  not  yet  clear  how  much 
this  uncertainty  affects  the  final  calculated  value  of  J , but  this,  and  poss- 
ible  methods  of  circumventing  the  problem,  will  be  the  subject  of  future 
investigation. 

20.  Both  the  proposed  method  of  J determination,  and  the  Wells  formula  for 

c 

5 , although  quite  amenable  to  calculation  by  hand,  involve  a fair  amount  of 
0 

arithmetic,  which  can  be  tedious  when  a large  number  of  specimens  is  involved. 

A small  computer  program  is  being  prepared  to  calculate  all  relevant  fracture 
parameters  from  the  basic  input  data  specified  in  the  Appendix. 
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EFFECT  OF  SLOW  CRACK  GROWTH  ON  THE  CORRELATION  BETWEEN  J & ABSORBED  ENERGY. 


Poi»4t  CmsplaccmEnIt 


Gau#e  Displacement  Curve  of  a 3pt  Bend  Specimen  0.25  < a/w  < 0.60 


From  the  specimen  measure: 

W,B,S  - overall  specimen  dimensions 

a^  - initial  crack  length  (to  end  of  fatigue  crack) 
a^  - final  crack  length  (to  end  of  slow  tearing) 

Draw  a line  0A  through  the  elastic  part  of  the  load/clip  gauge  deflection 

(Q/vJ  curve  (figure  5).  Draw  a horizontal  GF  at  load  Q„.  Record: 
e>  I 

5p(0G)  - final  load  (may  he  after  maximum  load) 

Q^(OH)  - limit  (net  section  yielding)  load 

V (OB)  - elastic  clip  gauge  displacement  (assuming  no  slow  crack  growth) 

V (OE)  - total  clip  gauge  displacement 
St 

Calculate: 


where 


TM  V 2 


2.9  - 4.6  (-£)  + 21.8  (^)2  - 37.6  (^)3  + 3b. 7 (-f ) 4 [A3) 


Calculate: 


2Ql(W  - af)  w 

w-  a^2  [af  + r(W  " af } + 


V . - j3  V 
gt  ge 


where  z is  the  height  of  the  clip  gauge  above  the  specimen  surface 
af. 

h(T^ 

e = “sr 

h(T) 

h(J)  = - 43.0  + 403-0  | - 107  3.2  (|)2  + 1162.8  (§)3 

r = 0.45  » a/W  <0.45 
0.4  > a/W  > 0.45 

The  value  of  J at  any  given  clip  gauge  displacement  is 


■QMPPHPU^M 


Notes 

a.  It  is  necessary  to  decide  whether  to  use  the  plane  strain  or  plane 

stress  value  of  S'  in  [Al]«  A possible  basis  for  cnoice  might  be  the 

criterion  for  determining  J,  validity  suggested  in  (26) 

X c 


b. 


c. 


d. 


Boy  &Y 

plane  stress  plane  strain 


[A9] 


[A4]  i3  not  applicable  if  significant  slow  era c.:  growth  occurs  before 
net  section  yielding.  can  be  calculated  before  net  section  yieldir 

in  the  absence  of  slow  crack  growth  by  substituting  ^ for  ^ in  [A4], 


[A6]  is  a polynomial  fit,  over  the  range,  0.25  < a/';T  < 0.6,  to  elastic 
compliance  data  given  in  (34).  Its  use  outside  this  range  is  not 
recommaided.  A slight  error  might  be  introduced  by  applying  this  to 
the  clip  gauge  as  opposed  to  the  mouth  opening  displacement,  but  this 
is  likely  to  be  small  since  (3  is  a ratio  not  an  absolute  quantity. 

The  values  of  r given  in  [A7]  are  based  on  finite  element  analysis. 
Alternative  values  may  be  used  if  experimental  evidence  dictates. 
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